Dynamical algebra and Dirac quantum modes 
in Taub-NUT background 

Ion I. Cotaescu * 

West University of Timi§oara, 
V. Pdrvan Ave. 4, RO-1900 Timi§oara, Romania 

Mihai Visinescu ^ 

Department of Theoretical Physics, 
National Institute for Physics and Nuclear Engineering, 
P.O.Box M.G. -6, Magurele, Bucharest, Romania 

February 1, 2008 



Abstract 

The 50(4, 1) gauge- invariant theory of the Dirac fermions in the 
external field of the Kaluza-Klein monopole is investigated. It is shown 
that the discrete quantum modes are governed by reducible represen- 
tations of the o(4) dynamical algebra generated by the components of 
the angular momentum operator and those of the Runge-Lenz opera- 
tor of the Dirac theory in Taub-NUT background. The consequence 
is that there exist central and axial discrete modes whose spinors have 
no separated variables. 

Pacs 04.62. +v 



1 Introduction 

In the relativistic quantum mechanics in the context of general relativity an 
important problem is how could be found the conserved operators which 
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should be involved in the sets of commuting observables defining quan- 
tum modes. From this point of view examples of geometries giving rise to 
large sets of conserved quantities are useful for understanding the connection 
among symmetries and conservation laws in general relativity. 

One of the most interesting geometries is that of the Euclidean Taub-NUT 
space since it has not only the usual isometries but also admits a hidden sym- 
metry of the Kepler type if a cyclic variable is gotten rid of 0, 0] . Moreover 
in the Taub-NUT geometry there are four Killing- Yano tensors among them 
the first three are covariantly constant and realize the quaternionic algebra 
related to the hyper-Kahler structure of the Taub-NUT manifold. The last 
one, that exists by virtue of the metric being of type has a non-vanishing 
field strength. These Killing- Yano tensors represent a certain square root 
of three Stackel-Killing tensors connected with the components of an ana- 
logue to the Runge-Lenz vector of the Kepler type problem since they give 
rise to constants of the geodesic motion (i.e. prime integrals) quadratic in 
4- velocities |, 

In other respects, the Taub-NUT metric is involved in many modern 
studies in physics. For example the Kaluza-Klein monopole was obtained 
by embedding the Taub-NUT gravitational instanton into five-dimensional 
Kaluza-Klein theory. On the other hand, in the long-distance limit, ne- 
glecting radiation, the relative motion of two monopoles is described by the 
geodesies of this space 0. The Dirac equation in the Kaluza-Klein monopole 
field was studied in the mid eighties j^. An attempt to take into account 
the Runge-Lenz vector of this problem was done in 0. We have continued 
this study showing that the Dirac equation is analytically solvable and 
determining the energy eigenspinors of the central modes. Moreover we de- 
rived all the conserved observables of this theory, including those associated 
with the hidden symmetries of the Taub-NUT geometry. Thus we obtained 
the Runge-Lenz vector-operator of the Dirac theory, pointing out its specific 
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In the present paper we investigate the consequences of the existence 
of the Runge-Lenz operator in the theory of the Dirac field in Taub-NUT 
background. We show that the dynamical algebras corresponding to different 
spectral domains are the same as in the scalar case but the representations of 
these algebras are different. Thus for each discrete energy level we obtain two 
irreducible representations of the o(4) algebra determining different quantum 
modes. This is a new phenomenon due to the spin terms of the Runge-Lenz 
operator since in the scalar case for each energy level we have only one 
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irreducible representation of the o(4) algebra 0. This conjecture offers us 
the opportunity to define new Dirac modes, called natural central or axial 
modes, the energy eigenspinors of which can be easily written using only 
algebraic properties. 

In Sec. 2 we briefly present the Killing vectors of the isometrics of the 
Taub-NUT background and the Killing tensors corresponding to the hidden 
symmetries of this geometry, pointing out that these are related to the four 
specific Killing- Yano tensors. The general form of the conserved operators 
of the Dirac theory is discussed in the next section while in Sec. 4 we write 
down the conserved operators arising from isometrics or hidden symmetries, 
presenting their algebraic properties. Sec. 5 is devoted to the study of the 
representations of the dynamical algebra of the discrete energy levels. These 
can be completely investigated since the action of the Dirac operators can be 
calculated at the level of two-component Pauli spinors where the structure 
of the corresponding operators is simpler. This method helps us to define 
the natural discrete Dirac modes in the next section where, in addition, we 
completely determine the form of the energy eigenspinors of these modes in 
terms of the spinors of the simple modes defined in Appendix. We work in 
natural units with h = c = 1. 



2 The Taub-NUT geometry 

The background of the gauge-invariant five- dimensional theory of the Dirac 
fermions in the external field of the Kaluza-Klein monopole 0] is the space 
of the Taub-NUT gravitational instanton with the time coordinate trivially 
added. Herein it is convenient to consider the static chart of Cartesian coor- 
dinates x'^, {fj,,!/, ... = 0, 1, 2, 3, 5), with the line element 

ds^ = gf.^dx^dx" = dt^ - ^df - V{dx^ + Aidx'f , (1) 

where dt^ = {dxY = {dx^Y + {dx'^Y + (dx^)^ is the usual Euchdean three- 
dimensional line element. This involves the Cartesian physical space coor- 
dinates = 1,2,3) which cover the domain D. The other coordi- 
nates are the time, x^ = t, and the Cartesian Kaluza-Klein extra-coordinate, 
x^ G -D5. Another chart suitable for applications is that of spherical co- 
ordinates, x' = {t, r, 6, 0, x), among them the first four are the time and 
the common spherical coordinates associated with x*, {i = 1,2,3), while 
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X + = —fix^. The real number /i is the main parameter of the theory 
which enters in the expression of the function 1/V{r) = 1 + /i/r. The unique 
non-vanishing component of the vector potential in spherical coordinates is 
= — cos 6'). In Cartesian coordinates we have 

divA = , 5 = rot 1 = u4 ■ (2) 

We note that the space domain D is defined such that 1/V > 0. 

The spacetime defined by (|I|) has the symmetry given by the isometry 
group Gs = S0{3) ® f^(l) <S> Tt{l) of rotations of the Cartesian space coor- 
dinates and and t translations. The ^75(1) symmetry is important since 
this eliminates the so called NUT singularity if x^ has the period Anfi. On 
the other hand, the 5*0(3) isometry transformations involve all of the space 
coordinates, and x^, since these are given by a non- linear representation 
of the 5*0(3) group. For this reason the corresponding Killing vectors, 
(z = 1, 2, 3), have the components /c^^-, = SijkX^ and /c^^^ = —eijkX^Ak—fix'^/r 
while the other two Killing vectors of the translation groups have usual con- 
stant components, k'^Q^ = and k'^^^ = 6^. 

The Taub-NUT geometry possesses four Killing- Yano tensors, f^^^ {i = 
1, 2, 3) and f^, of valence 2 (which satisfy ff.^ = -fufi, and f^u;a + f,i<7;u = 0). 
The first three, 

fi^) ^ /\ ^ A + e,,ke' A e'^ , (3) 

are rather special since they are covariantly constant (with vanishing field 
strength). The fourth Killing- Yano tensor of the Taub-NUT space, 

/^ = --/« + ^^.,.e^■Ae^ (4) 
r /iV/ 

is not covariantly constant having the non-vanishing components of its field 
strength 

2r2 

/r,;, = ^sin^, (5) 

which are completely antisymmetric in r, 9 and (f). 

The hidden symmetries of the Taub-NUT geometry are encapsulated in 
the non-trivial Stackel-Killing tensors k(^i)^^,{i = 1,2,3). They can be ex- 
pressed as symmetrized products of Killing- Yano tensors (|^) and (^) [^]: 

ki^),u = -^{f\xf%' + flxf%') + + %).%),)• (6) 
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In fact only the product of Killing- Yano tensors f^^^ and leads to non- 
trivial Stackel-Killing tensors, the last term in the r.h.s. of (j^) being a simple 
product of Killing vectors. This term is usually added to write the Runge- 
Lenz vector in the standard form of the scalar (classical, Schrodinger or 
Klein-Gordon) theory 0, |, IH. 



3 The operators of the Dirac theory 

The theory of the Dirac field in Cartesian charts of the Taub-NUT ge- 
ometry takes the simplest form if one considers the local frames given by 



pentad fields, e{x) and e(a;), as defined in |0. Their components, have 
the usual ort honor malizat ion properties and give the components of the 
metric tensor, g^i, = V&p^'^^u g^'^ = rj°'^e^e^i^. In our notation 0, 
7] = diag(l, —1,-1,— 1,-1) is the fiat metric which raises or lowers the 
hated indices (ranging from to 5). Its gauge group, G{ri) = 50(4,1), 
has as universal covering group, G{ri), a subgroup of the group SU{2,2) 
carried by the space of four-dimensional Dirac spinors. Therefore, the five 
matrices 7", that must satisfy {7", 7^} = 27]°'^, can be defined in terms of 
the standard Dirac matrices |]T^, |T3[ as 7° = 7°, 7* = 7*, {i = 1, 2, 3) and 
7^ = 27^. These are self-adjoint with respect to the usual Dirac conjugation, 
i.e. 7" = 70(7")t^o = and give the covariant basis generators of the 
group G{ri), denoted by S""^ = ^[7°, 7'^]/4. 

The Dirac field ip of mass M, defined on the space domain D x D^, has 
the gauge-invariant action 0, |^ , 



J (fx^ 1^ [^7^V^^ - (V^V^)7^V^] - MV-^} (7) 



where Vq are the components of the spin covariant derivatives with local 
indices P, 



iVVPi + '-V^/Vsi.unBk , V5 = ^7^5 - l^V^/vt* -B. (8) 



'7v"'~2 

These depend on the momentum operators Pi = —i{di — Aid^) and P5 = 
—185 , which in Taub-NUT geometry obey the commutation rules [Pi, Pj] = 
iSijkBkP^ and [Pi, P5] = 0. The spin matrices which give the spin connection 
are defined by 

S: = 5,, + i7V, S, = lei,,S^'. (9) 
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The action (0) leads to the Dirac equation, Dip = Mip, given by the Dirac 
operator ^ 

= z7°9t-v^7-P-^75p5-^Vv^75S*-5, (10) 

yV 2 



having the usual time-dependent term |]T5| and a static part, T>s. 

In the standard representation of the Dirac matrices (with diagonal 7° 
||15||) the Hamiltonian operator is p[ if = if + 7°M. The massless Hamilto- 



V-K* ^ 



man, 

H = 7°P. = 1 1_ "VV\, (11) 
VVtc 

is expressed in terms of the operators vr = ap — iV'^P^ and tt* = ap+iV^^P^ 
where ap = a ■ P involves the Pauli matrices, cTj. These Pauli operators give 
the space part of the massless Klein-Gordon operator as, 

A = Vn*n = VP\ yP^^ . (12) 

In what follows, we are interested especially by the form and the action 
of the conserved operators of the Dirac theory which, by definition, are the 
operators that commute with the Hamiltonian operator H. Since the mass 
term is simple and does not rise difficulties to keep the study of the operators 
as simple as possible we leave it aside. We shall consider the massless case 
(M = 0) and the Hamiltonian operator reduces to ([TT|) without affecting the 
symmetries and conserved observables of the Dirac theory. The conserved 
operators depend on Pauli operators obeying several conditions requested by 
the above definition. Two types of Dirac operators are important for our 
further developments. The first one are the Q-operators introduced in M as 



Q(X) = {H, 





(13) 



where Qo = iVg = ij'^H. The linear mapping Q associates Dirac operators 



to the Pauh operators, vr, vr*, ap, ai = a ■ L, ar = a ■ x/r, etc. ||T6[. The 
remarkable property of the Q-operators is that if [X, A] = then we can 
write 

QiX) Q{Y) = HQ{XY) , Q{Y) Q{X) = QiYX)H , (14) 
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for any other Pauli operator Y . Hereby it results that Q{X) commutes with 
H since H = Q{1). Moreover, one can verify that if both the operators, X 
and Y, commute with A and commute or anticommute between themselves 
then Q{X) and Q{Y) commute or anticommute each other. An important 
property is that all of the Q-operators anticommute with 7°. Another type 
of conserved operators are of the diagonal form T = diag (T*^^\ T*^^)). Then 
the condition [T, H] = requires the Pauli operators T^^^ and T^^^ to satisfy 



m 



T^^^Wn = WnT^^^ , Vtt*^T^^'^ = T^^^Vn*^ . (15) 



Hereby it results [T^^\ A] = and the useful relations 

TQ{X) = Q(rWx) , Q{X)T = Q(Xr«) , (16) 
that can be used for deriving new algebraic properties. 



4 Symmetries and conserved observables 

A class of conserved operators are the generators of the operator-valued rep- 
resentations of the isometry group carried by the spaces of physical fields. 
In the scalar case, these are orbital operators defined up to the factor —i as 
the Killing vectors fields associated to the isometrics 0]. For the fields with 



spin the conserved generators get, in addition, specific spin terms determined 
by the form of the Killing vectors and the gauge fixing. Thus, in the case of 
the Dirac field (in four or five dimensions), for each Killing vector k one can 



write the operator 17, 18 



X, = -ik^elV^ + \k^;ue%elS''^ (17) 

which commutes with the Dirac operator T). 

In Cartesian coordinates and our gauge fixing, the f/t(l) generator is 
—idti while the U^{1) one is P5. In spherical coordinates it is convenient to 
replace P5 with Q = — yuPs = —idx- The other three Killing vectors give the 
5*0 (3) generators which are the components of the whole angular momentum 
J' = L + S as in the flat spacetimes. The difference is that here the orbital 
angular momentum, 

L = XX P- U-P5 , (18) 
r 
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depends on P5 since the S0{3) isometries are non- linear transformations. 
The consequence is that the irreducible representations of the o(3) algebra 
generated by Lj are similar with the linear ones but with a supplementary 
restriction upon the angular quantum numbers [0. On the other hand, the 
operators J7i = diag(Jj, Jj) have diagonal form where Jj = L, + ai/2 are 
just the angular momentum operators of the Pauli theory. For this reason 
J'i commute with T> and 7'^ and, therefore, they are conserved, commuting 
with H. Moreover, they satisfy the canonical commutation rules among 
themselves and with the components of all the other vector operators (e.g. 
coordinates, momenta, spin, etc.). 

However, there are other types of conserved operators directly related to 
the specific geometric objects of the Taub-NUT geometry, as Stackel-Killing 
and Killing- Yano tensors. The Killing- Yano tensors give rise to conserved 
observables defined as Dirac-type operators of the form 

Qf = fo^pT^^ - \nx)Y{x)^\x)f,^,,{x) (19) 

where f^^^ = ffj.ue-%e^i^ and 7'^(x) = e^-^"-. According to an important result of 
Ref. ||18||, these operators anticommute with the Dirac operator V. 



Starting with the first three Killing- Yano tensors, from Eq.(|T9D, after 
some algebra, we obtain the Dirac-type operators 

= = 2(^0 ' (20) 

which anticommute with Qq and 7*^, commute with H and obey 

QiQj = SijH'^ + iSijkQkH , [J7i, Qj] = iSijkQk ■ (21) 

Hereby we find the = 4 super algebra 

{Qa,Qb} = '25abH\ a, 5,... =0,1,2,3. (22) 

The corresponding Dirac-type operator of the fourth Killing- Yano tensor, 
calculated with the general rule ([l9|) , is 



g^ = -Q(a.) + ^ \ : (23) 




where X = a ■ {x x P) + 1 = + 1 + ^(TrP^ is the operator studied in 
||19||. Other equivalent forms are given in |jTO|. Again one can verify that 



commutes with H and anticommutes with Qq and 7 
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Of a special interest is the Runge-Lenz operator of the Dirac theory as- 
sociated to the Kilhng tensor k'^'^. This can be constructed with the help of 
the conserved Dirac-type operators generated by the Killing- Yano tensors. 



One defines first the vector operator of components |T0 



(24) 



which commutes with H since the operators and Qi are commuting with 
H. Consequently, N has diagonal form, its diagonal blocks, N^^^ and N^'^\ 
obeying the general conditions (15). One finds that the first block which 
commutes with A, 



contains not only the orbital Runge-Lenz operator. 



/i X 



X 



K = -{P X L-Lx P)- --A + /i-P; 



2 r 



5 ; 



(25) 



(26) 



but a spin term too. Furthermore, by taking into account that the compo- 
nents of K commute with A and satisfy Q 



A 



(27) 
(28) 



one can calculate the algebraic properties given in |I0[ among them the com- 
mutation relations 



[iVj, Nj] — iSijkJk^'^ + 2^ijkQkH , T'^ — Ps^ 



W , (29) 

suggested us to define the components of the Runge-Lenz operator of the 



Dirac theory, /C, as follows |T0 



1 



Then one obtains the commutation relations 



[/C„ if] = 0, 
[/C,, P5] = 0, 



[A^j, 3^ j\ i^ijk^k ; 
[A^j, Qj] i^ijkQk^ : 



(30) 



(31) 
(32) 
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and 

[/Cj, K,j\ = iSijkJk^ ■ (33) 

Since there are no zero modes the operator H is invertible such that our 
definition of the Runge-Lenz operator is correct.. It is worthy to note that /Cj 
are diagonal commuting with 7° which means that they are also conserved 
in the massive case ([/C,, iJ] = 0). 



5 Dynamical algebra 

The large collection of conserved observables we have obtained will help us 
to select many different complete sets of commuting observable which should 
define static quantum modes, with a given energy > 0. In other respects, 
since P5 commutes with all conserved observables, its eigenvalue, g, will play 
the role of a general parameter. When one uses the operator Q instead of P5 
then we have to consider as parameter its eigenvalue q = —^q. 

In these conditions we can re-scale the Runge-Lenz operator in order to 
recover the dynamical algebras o(4), o(3, 1) or e(3), corresponding to different 
spectral domains of the Kepler-type problems p|. If we define 

{T^^Ki for ;U < and E < \q\ 
/Cj for any /i and E = \q\ (34) 

iJ-'^^JCi for any /i and E > \q\ 

then the operators J'i and TZi {i = 1,2,3) will generate either a representation 
of the 0(4) algebra for the discrete energy spectrum in the domain E < \q\ or 
a representation of the o(3, 1) algebra for continuous spectrum in the domain 
E > \q\. The dynamical algebra e(3) corresponds only to the ground energy 
of the continuous spectrum, E = \q\. 

The manipulation of the dynamical algebras seems to be difficult in view 
of the complicate structure of their generators but in fact the calculations can 
be done at the level of Pauli operators because of the special properties of the 
common eigenspinors of H and P5. These have the form UE,q = {u^^q, u^ej)'^ 
depending on two-component Pauli spinors which satisfy the equations @] 

A«S = (35) 

= E-'Wnu^^l, (36) 
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equivalent with the eigenvalue problem {H — E)UE,q = 0. This means that 
u^Eq may be any solution of the scalar equation (|35|) which is just the static 
Klein-Gordon equation in Taub-NUT geometry. Moreover, it is clear that 
u^Eq completely determined by Eq. (|36D if the form of u'Eq known. Since 
the Klein-Gordon equation is scalar, the form of the spinor u~Eq s-nd, there- 
fore, that of UE,q is strongly dependent on the choice of the other observables 
included in the complete set of commuting operators which defines the quan- 
tum modes. Technically speaking, we have to start with the first Pauli spinor 
of the form Wg^*- ~ X^o- /o-lo- where ^o- are the usual Pauli eigenspinors of po- 
larizations a = ±1/2 while fa are scalar solutions of Eq. ([35|) . The next step 
is to precise this linear combination with the help of the other conserved 

(2) 

operators of the complete set of commuting operators and to write 
cording to (pG)) . The action of these operators can be easily calculated using 
the mentioned general properties (ITB) of the diagonal operators or the struc- 



ture of the off-diagonal ones P]. For example, it is not difficult to show that 
the action of the Runge-Lenz operator upon the energy eigenspinors is 

lCUEq=i ^ ^ , K' = K + F-, (37) 

\^ E-'VVnK'u^E^q J 2 ' ^ ^ 

which means that this reduces to the action of the new operator K' upon the 
first Pauli spinor. We say that K' is the Pauli operator associated with K,. In 
fact, because of the central role played here by u"E,qi this property is general, 
each conserved observables having its own associated Pauli operator. Thus 
any problem of the Dirac theory in Taub-NUT background can be rewritten 
in terms of Pauli operators acting upon the first two-component spinors of 
the spinors UE,q- 

In order to illustrate how works this mechanism, let us study the rep- 
resentations of the o(4) dynamical algebra of the discrete quantum modes 
of the Dirac field. As mentioned, these arise in the domain E < \q\, only 
when fi < 0, and have the same energy levels as the scalar Klein-Gordon (or 
Schrodinger) equation. 



2 2 



(38) 



These are labeled only by the principal quantum number n which takes all 
the integer values larger than |g| 0. The quantization rule of the scalar 
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modes is related to the irreducible unitary finite-dimensional representations 
of the o(4) dynamical algebra generated by Li and the components of the re- 
scaled Runge-Lenz vector, i?, = F~^Ki. Since the second Casimir operator 
of this algebra is C2 = L ■ -R = 0, the first one, Ci = + B?, can take only 
the eigenvalues ci = — 1 (with n > \q\) which give the energy levels (^8|). 
In terms of su{2) weights these representations of o(4) ~ su{2) x su{2) are 
denoted by (^,^). 

In the Dirac case the dynamical algebra is the same but its representa- 
tions are generated by J'i and TZi = JF~^/Ci. According to Eq. (|37[) , these 



representations are equivalent with those generated by the associated Pauli 
operators Jj and R'^ = F~^K[ acting upon the first Pauli spinor. However, 
since Jj = Li + ai/2 and R'i = Ri + ai/2, we draw the conclusion that the 
Dirac discrete modes are governed by the reducible representation 

n — 1 n — 1\ (\ \ (n n — \\ fn n — 1\ , , 

— — )s(5.o) = (^— )e(2-i— ). (39) 

The Casimir operators, C[ = + R' and €'2 = J ■ R' take now the eigenval- 
ues, c'l and C2, such that c'l — = — 1 for both irreducible representations 
while C2 = (2n + l)/4 for the representation (|, ^) and 4 = -(2?t, - l)/4 
for the representation (| — 1, ^^). This suggests us to use the new Pauli 
operator 

C = 2C^ - 1/2 = aR + (jL + l, (jR = a-R (40) 

in order to distinguish between the irreducible representations resulted from 
the decomposition (^). The advantage is that this has the simplest eigen- 
values, c = ±n. 



6 Discrete quantum modes 

Since the energy levels En are degenerated, we need to use complete sets 
of commuting operators for determining quantum modes. Fortunately the 
set of conserved observables is large enough to offer us many possibilities 
of choice. An appropriate option is to consider the natural modes involving 
only one irreducible representation of (PSf). The complete sets of commuting 
observables of these modes must include the operators H, P5 and the new 
operator C = 2J -71—1/2 associated with the Pauli operator C defined above. 
Then only two more operators we need for defining these Dirac quantum 
modes. 



12 



We say that the set {H, P5, C, JT"^, J3} defines the natural central modes. 
Its common eigenspinors, Un,q,c,j,m , correspond to the eigenvalues E^, q, c, 
j{j + 1) and rrij. Another possibihty is to take the set {if, P5, C, 71^, jTs} 
of the natural axial modes the common eigenspinors of which, Un,q,c,mr,mj, 
correspond to the eigenvalues En, q, c, rrir and mj. On the other hand, 
our previous results O, O lead to the conclusion that there are simple 



modes having eigenspinors with separated variables. These are central and 
axial simple modes defined by the sets {H, P5, JT"^, jTa, Q(ctl + 1)} and 
{H, P5, 7I3, jTs, Qs} respectively, as it is shown in Appendix. Since neither 
Q((7l + 1) nor Q3 do not commute with C, it results that the natural modes 
do not have eigenspinors with separated variables. Therefore these must be 
linear combinations of the spinors of simple modes. Then it is interesting to 
try to write down these linear combinations using only the algebraic method 
based on the relation between Dirac and Pauli operators. 

According to this method, we observe the first Pauli spinors, u^^^^^j^rnj^ of 
the eigenspinors Un,q,c,j,mj of the natural central modes must be the eigen- 
spinors of the set {A, P5, C, J^, J3} corresponding to the eigenvalues E^, q, 
c, j{j + 1) and rrij. On the other hand, the superalgebra 



{an, aL + l} = 0, {anY = C, - P 



(41) 



allows us to demonstrate that the first Pauli spinors of the eigenspinors of 
the simple central modes ( [A.l| ) have the remarkable property 



(i)± 

^n,q,j, rrij 



nl/2 



u. 



{1)T 
n,q,j, rrij 



(42) 



Hereby it results that the Dirac eigenspinors of the natural central modes 
can be expressed as 



Un,q,c=±n,j ,% 



Mh±{j + VU^,q,J,m, 



+ \lnT [J + I) U^,q,j,m, 



(43) 



In the same way the eigenspinors of the natural axial modes can be 
expressed as linear combinations of the spinors of the simple axial modes 
presented in Appendix. As in previous case the calculations reduce to the 
eigenspinors of the set of Pauli operators {A, P5, C, P3, J3} associated to 
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that of the Dirac operators of the natural axial modes. Using the identity 
{C, 0-3} = 2{R'^ + J3) it is not difficult to show that the eigenspinors of the 
natural axial modes are 



m,r,mj 



'In 



±\Jn ± \mr + mjl U^^^ 



mr,mj 



(44) 



Thus we see that in the Dirac theory the algebraic method offers us 
the mechanisms of constructing new quantum modes having no separated 
variables. We can say that this method and that of separation of variables 
complete each other, helping us to find many types of different quantum 
modes related among themselves. 



7 Concluding remarks 

The aim of this paper was to complete the previous studies of the Dirac 
equation in Taub-NUT space p, [1^. In the case of the Dirac equation the 
existence of the Killing- Yano tensors in the Taub-NUT geometry allowed us 
to construct some Dirac-type operators the anticommutators of which gener- 
ate operators that commute with the standard Dirac one. We are attempting 
to find non-trivial operators of the Dirac theory connected with the hidden 
symmetry of the Taub-NUT space. These operators are constructed in Sec. 4 
and they represent the quantum analogue of the classical Runge-Lenz vector 
from the Kepler problem. Thus we get an example of gauge-invariant theory 
of Dirac fermions in a geometry with high manifest or hidden symmetries 
for which we can write all the corresponding conserved operators. The main 
consequence is that we can explicitly use the dynamical algebra for defining 
new quantum modes. 

In the literature there is a detailed description of the classical geodesic mo- 
tion for scalar particles and a quantum treatment through the Schrodinger 
equation in Taub-NUT space [jl], 0, |12|. On the other hand, the pseudo- 



classical limit of the Dirac theory of a spin— | fermion in curved spacetime 
is described by supersymmetric extension of the ordinary relativistic point 
particle 12^]. In the pseudo-classical models the spin degrees of freedom are 



characterized in terms of anticommuting Grassmann variables. The con- 
stants of motion related to the symmetries of the manifold generally contain 
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spin-dependent parts It is interesting that these specific spin contribu- 
tions are similar to the spin terms of the conserved operators of the quantum 
theory. A notable example is the spin contribution to the Runge-Lenz vec- 
tor which re-confirms the result from pseudo-classical approach [Q. More- 
over, since both the pseudo-classical and the quantum theory are completely 
solvable, we have the opportunity to compare their physical meaning in all 
details. 



Appendix A: The simple discrete modes 

The energy eigenspinors of the simple central modes have the form 

TT± — n,q,j,mj J n,qj ^ q,j,mj \ (AT\ 

\ "'n,q,j,mj ] \ q,J,mj ~ yn,q,j ^ q,],mj J 

where the radial functions Z"^, and depend only on r while are the 
two-component spherical spinors that solve the common eigenvalue problems 
of the Pauli operators Q = — /uPs, J^, J3 and + 1 for the eigenvalues 
q = — yug, j{j + 1), rrij, and ±(j + |), respectively 0. We observe that 
only the first Pauli spinor has separated variables. Its radial function, /, 
is a solution of the radial Klein-Gordon equation while the other two radial 
function arise from (|36|). The quantum numbers of these modes must satisfy 
the selection rule lol < ? ' + ^ < n. Thus we find that u^^^, ^ is the common 
eigenspinor of the set {A, P5, J^, J3, (7^+1} corresponding to the eigenvalues 
-^n' jU + + |) respectively. This means that (|A.1| ) is 

the common eigenspinor of the operators {H, P5, jf^, jTs, Q{aL + 1)} whose 
eigenvalues are En, q, j{j + 1), mj and ±_E„(j + |). 

The simple axial modes can be constructed in a similar way. We start 
with scalar common eigenfunctions fn,q,m,m of the commuting operators A, 
P5, P3 and L3, corresponding to the eigenvalues E^, q, m and m. These 
eigenfunction can be easily calculated in parabolic coordinates as the axial 
solutions of the Schrodinger equation , with the unique difference that the 
quantity 2E from the Schrodinger case must be replaced here by E"^. Fur- 
thermore, we construct the eigenspinors of the simple axial modes, Un^^^nir,mj: 
by choosing their first Pauli spinors as 

=±1/2 (A.2) 
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and calculating from ( P^ ) the second Pauli spinors that have no more sepa- 
rated variables. It is clear that ( |A.2|) is the common eigenspinor of the set 



of Pauli operators {A, P5, i?3, J3, 0-3} associated with the Dirac ones from 
the set defining the simple axial modes. The eigenvalues of these operators 
are E^, q, mr = m + a, rrij = m + a and a. Consequently, U^^^^^^^^. are the 
common eigenspinors of the set {H, P5, 7^3, jTs, Qs} corresponding to the 
eigenvalues g, m^, rrij and EnCr. 
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